Magnetic-field-dependent quasiparticle energy relaxation in mesoscopic wires 
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In order to find out if magnetic impurities can mediate interactions between quasiparticles in 
metals, we have measured the effect of a magnetic field B on the energy distribution function f(E) 
of quasiparticles in two silver wires driven out-of-equilibrium by a bias voltage U. In a sample 
showing sharp distributions at B = 0, no magnetic field effect is found, whereas in the other sample, 
rounded distributions at low magnetic field get sharper as B is increased, with a characteristic field 
proportional to U. Comparison is made with recent calculations of the effect of magnetic-impurities- 
mediated interactions taking into account Kondo physics. 

PACS numbers: 73.23.-b, 72.15.Qm, 72.10.-d, 71.10.Ay 
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The understanding of the phenomena which, at low 
temperature, limit the extent of quantum coherence in 
electronic transport and allow the quasiparticles to ex- 
change energy is presently an important issue in meso- 
scopic physics. There is indeed a discrepancy between 
the theory 0, which predicts that Coulomb interactions 
provide the dominant mechanism for decoherence and 
for energy exchange, and measurements of the coherence 
time 0, @ or of energy exchange rates 0, 0, IE m 
numerous metallic samples. This discrepancy has been 
attributed either to a flaw in the theory or to the 
presence in these samples of other mechanisms involv- 
ing the scattering of electrons by undetected two-level 
systems or magnetic impurities. It has been indeed re- 
cently predicted that even a minute concentration of such 
scatterers would result in sizeable energy exchange if the 
Kondo effect occurs [E IE ■ Whereas the limitation of 
quantum coherence by the Kondo effect is widely known 
, its efficiency for mediating energy exchange between 
quasiparticles had not been anticipated. In the case of 
magnetic impurities, a significant weakening of this effec- 
tive electron-electron interaction is furthermore predicted 
when a large magnetic field is applied ^E|- I n order to 
test these new predictions and more generally to under- 
stand inelastic processes in mesoscopic conductors, we 
have investigated the magnetic field dependence of the 
energy exchange rate in mesoscopic wires. 

The samples are wires connected to reservoirs biased 
at potentials and U (see Fig. 1). The energy distribu- 
tion function in the middle of the wire, f{E), depends 
on the ratio of the typical interaction time Ti nt and the 
diffusion time of quasiparticles td = L 2 /D. If Tint 3> t~d, 
interactions can be neglected and f(E) is the average of 
the Fermi functions in both reservoirs, which have elec- 
trochemical potentials shifted by eU. In the experimental 
situation where fcgT <C eU, f(E) is then a two-step func- 
tion. In the opposite limit Tint <c to , local equilibrium is 
achieved at each coordinate along the wire, and f(E) is 
a Fermi function at a temperature given by the balance 
between Joule heating and electronic heat conductivity 
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FIG. 1: Top: layout of the experiment: a wire is connected to 
two large electrodes biased at potentials and U. A resistive 
probe electrode (in grey) forms a tunnel junction with the 
wire. At equilibrium (U = 0), the differential conductance 
dI/dV(V) of this junction displays a dip at zero bias, due to 
Coulomb blockade of tunneling (left). When the quasiparti- 
cles of the wire are driven out-of-equilibrium by a finite volt- 
age U (right), their energy distribution function f(E) depends 
on the interaction rate between quasiparticles. In the absence 
of interactions, f(E) is a two-step function and dI/dV(V) 
presents two dips (solid lines). With strong interactions, f(E) 
is rounded, and dl /dV(V) presents a broad dip (dashed lines). 



to the reservoirs: this is the "hot-electron" regime [l3j. 
The intermediate regime is of interest for experiments be- 
cause the precise shape of f(E) and its dependence on U 
are characteristic of the interaction rate and of its energy 
dependence |4J. 

At zero magnetic field, the distribution function 
f(E) can be inferred from the differential conductance 
dl / dV (V) of a tunnel junction between the central part 
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of the wire and a superconducting (aluminum) probe 
electrode biased at potential V Q. In magnetic fields 
larger than the critical field B c ~ 0.1 T of the supercon- 
ducting electrode, another method is required. Here, we 
have taken advantage of the nonlinearity of the current- 
voltage characteristic of a tunnel junction placed in series 
with a resistance R. When both electrodes of the junc- 
tion are in the normal state and at thermal equilibrium, 
the differential conductance dI/dV(V) displays a dip at 
V = (see Fig. 1), due to the Coulomb blockade of tun- 
neling 14]. Assuming that the two electrodes have dif- 
ferent distribution functions / and / rc f, the differential 
conductance reads: 

§(Y) = ±JdEf(E) JdeP(e) (1) 

x ^((f Ict (E + eV + e) - f Icl (E + eV - e)) 

where Rt is the tunnel resistance of the junction, and 
P(e) = J e J(t)+iet/h p ro b a biiity for an electron 
to tunnel through the barrier while releasing to the envi- 
ronmcnt an energy e, J(t) = f ^ aR ff")l ^Zlj^ 
with Z(uj) = \/{l/R+jCw), C the junction capacitance, 
Rk = h/e 2 k, 25.8 kSl the resistance quantum and T the 
environment temperature. In the case where the distri- 
bution function f(E) presents two steps, as in Fig. 1, and 
/ rc f is a Fermi function at temperature T, one obtains, by 
linearity, two dips in dI/dV(V) at V = and V = — U. 
In contrast, in the hot electron regime, dI/dV(V) dis- 
plays a broad dip centered at V = —U/2 (see Fig. 1). 
In the experiments, a large series impedance at the rel- 
evant frequencies (up to about 50 GHz) was obtained 
by designing the probe electrode as a long, narrow and 
thin aluminum electrode (25/imx 150nmx 12nm), which 
presents a resistance R ~ 1.5 kS7 in the normal state. 

We present here the results obtained on two silver sam- 
ples in which the distribution functions found at B = 
were extremely different. The samples were obtained 
from nominally five- nines-purity (99.999%, sample #1) 
and six-nines-purity (99.9999%, sample #2) source ma- 
terial. For both wires, the length and cross-section area 
are L = 20 /itm, S = 100 nm x48 nm. The diffusion con- 
stants D = 196 and 215 cm 2 /s respectively, were deduced 
from the low temperature resistance. The tunnel resis- 
tances Rt (167 and 102 kSl) and the capacitances C (0.8 
and 0.9 fF) of the junctions, as well as the environment 
resistances R (1.34 and 1.65 kO), were obtained from fits 
with Eq. (JTJ of dI/dV(V) measured at B = 0.3 T and 
[7 = 0. We have checked that these curves do not change 
with B when B > B c . 

At low magnetic field and low temperature, the probe 
electrode is superconducting. Its impedance is purely 
imaginary at frequencies lower than 2A/h ^^]. It re- 
sults that for eV G [-3A + U, 3 A] Coulomb blockade 
only leads to a reduction of the differential conductance, 
which is multiplied by a factor exp(— ^ 2B gf ) 
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FIG. 2: a) Symbols: distribution functions f(E) at U = 
0.15 mV and zero magnetic field in samples #1 and #2, ob- 
tained by deconvolution of dI/dV(V) with the probe electrode 
in the superconducting state. Solid lines: fits with theory in- 
cluding the effect of Kondo impurities (see text), b) Dashed 
line: calculated Coulomb blockade signal dl /dV(V) using the 
measured f(E) at B — 0. Symbols: measured dI/dV(V) at 
U = 0.15 mV, with B = 0.3 and 1.2 T, the probe electrode 
being in the resistive state. Solid line: fits with theory. 



~ 0.9. Numerical deconvolution of dI/dV(V) is there- 
fore possible, and the distribution functions obtained at 
U = 0.15 mV are presented in the top of Fig. 2 for both 
samples. Whereas f(E) is close to a double-step function 
in sample #2, it is much more rounded in sample #1, in- 
dicating that the energy exchange rate is much larger 
in the latter, since the diffusion times are very similar 
(td = L 2 /D ~ 20 ns). In the bottom of Fig. 2, we plot 
the calculated R T dI /dV(V) using formula (Q) with f(E) 
the distribution function measured at B = (dashed 
curves), and present the measured curves for B = 0.3 T 
and B = 1.2 T (symbols) In sample #2, the mag- 
netic field has no visible effect. Note however that the 
distribution functions are so close to a double-step that 
the experiment is not sensitive enough to detect a pos- 
sible slight reduction of the energy exchange rate with 
B. In contrast, in sample #1, the rounded dip at zero 
field is replaced at 1.2 T by a double-dip, showing that 
the energy exchange rate has been reduced. Figure 3 
shows the evolution of dI/dV(V) with magnetic field, 
from 0.3 T to 1.5 T by steps of 0.3 T, for U = 0.1, 0.2 
and 0.3 mV. A similar behavior is observed at all values 
of U; the low-field broad conductance dip at B = 0.3 T 
tends to be replaced at large fields by a double-dip struc- 
ture. In particular, the crossover field at which dI/dV{V) 
is nearly constant over a broad voltage range is 0.6 T at 
U = 0.1 mV, 0.9 T at U = 0.15 mV (not shown), 1.2 T 
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FIG. 3: Magnetic field effect in sample #1: differential con- 
ductance dI/dV(V) at U = 0.1, 0.2 and 0.3 mV, for B ranging 
(from bottom to top) from 0.3 to 1.5 T by steps of 0.3 T. Suc- 
cessive curves have been vertically offset by steps of 0.033, for 
clarity. 

at U = 0.2 mV, and 1.5 T at U = 0.25 mV (not shown), 
hence presenting a linear increase with U. The compari- 
son of the raw data on sample #1 and sample #2 in Fig. 2 
already allows to conclude that sample #1 presents an 
extra interaction which can be strongly reduced by ap- 
plying a magnetic field. 

We now compare the experimental data with theoret- 
ical predictions. The distribution function is calculated 
by solving the stationary Boltzmann equation in the dif- 
fusive regime 0, 0] : 

^ 92f Q X x 2 E) = -4o U (^^{/})+^ou(^^{/}) (2) 

where X™ n and 1°"^ are the rates at which quasiparticles 
are scattered in and out of a state at energy E by inelastic 
processes. Assuming that the dominant inelastic process 
is a two-quasiparticle interaction which is local on the 
scale of variations of the distribution function, 

% = J dedE'K (e) f% +e (l - f%)f%,(l - /§,_ e ) (3) 

where the shorthand f% stands for f (x,E) . The out- 
collision term X™^ has a similar form. The kernel func- 
tion K (e) is proportional to the averaged squared in- 
teraction between two quasiparticles exchanging an en- 
ergy e. Coulomb interactions lead, in diffusive wires, to 

K(e) = k/s 3 / 2 Q, where k = {ir^tf^vpS)- 1 with 

vp the density of states at the Fermi level |l9j. The 
B = data for sample #2 can be well fit using this term 
with k = 0.12 ns _1 meV" 1 ^ 2 , of the same order of magni- 
tude as the theoretical value 0.07 ns _1 mcV -1 ^ 2 [13], and 
a term of lesser importance describing phonon emission 
|2lj . The B = data for sample #1 can be fit simi- 
larly, with k — 2.4 ns _1 meV ^ 2 , however the reduction 
of the energy exchange rate with B indicates that an 
extra process is present at B = 0. We have in the fol- 
lowing fixed k to the best fit value obtained for the large 
field, low U data, where the B— dependent interaction 



has essentially vanished: k = 0.5 ns _1 meV -1 ^ 2 [2(ij], The 
remaining part of the energy exchange rate was fit with 
the theory of Goppert et al. 0, H3| , which accounts for 
the effective interaction in the presence of a concentra- 
tion c of spin-i impurities, with Kondo temperature Tk, 
gyromagnetic factor <?, and coupling constant J between 
quasiparticles and magnetic impurities. The Kondo ef- 
fect is included in this calculation, under the assumption 
that the distribution functions are not too sharp, leading 
to a renormalization of J depending on the distribution 
function itself. The corresponding inelastic integral can 
also be written in the form of Eq. J2J), but with a K (e) 
function depending on the energies E and E' and on /. 
At zero magnetic field, the effect of this term on f(E) is 
similar to that of a phenomenological kernel K (e) oc 1/e 2 
as found in For compatibility with phase-coherence 
time measurements (see below), the Kondo temperature 
was fixed at Tk = 40 mK, which is the Kondo tempera- 
ture of Mn in Ag. As shown by solid lines in Fig. 2 and 
in Fig. 3, the data can be accuratel y re produced using 
c = 17 ppm, g = 2.9 and v. J =0.08 23] . Note however 
that according to material analysis of the silver batch 
used to fabricate sample #1, no magnetic impurity was 
present in the source at the level of 1 ppm. Since in some 
samples made out of the same batch the intensity of the 
energy exchange rate measured at B = was found to 
be up to 4 times smaller, pollution of the sample dur- 
ing fabrication might, at least partly, explain this large 
impurity concentration. 

The impurity concentration deduced from the fits of 
f(E) must be further compared with the one obtained 
from the analysis of measurements of the phase coher- 
ence time in long wires fabricated previously with the 
same source materials. We have extracted the phase 
coherence time t v from the magnetoresistance of wires 
much longer than the phase coherence length, using the 
weak localization theory. In samples made of 6N pu- 
rity Ag, t v (T) = A/T 2 / 3 from 1 K down to 40 mK, 
with A = 2.25 ns K 2 / 3 , in reasonable agreement with 
the theory of Coulomb interactions in disordered wires 
(A thoory = 3.00 ns K 2/3 ). At T = 40 mK, t v = 18 ns. 
In samples made of 5N silver, t v (T) does not vary be- 
tween T = 200 mK and 40 mK, where we find r v = 2 ns. 
This behavior can be attributed to the presence of mag- 
netic impurities, with concentration c, spin s and Kondo 
temperature Tk, which lead to a spin- flip rate described 
by El S3 7./(T) = (c/Tthv) n 2 s(s + l)/(7r 2 s(s + 
1) + \t?(T/Tk))- The resulting phase coherence time 
t v (T) = \/{T 2 / 3 /A + 27„/ (T)) shows very little varia- 
tion between 40 mK and 200 mK and describes precisely 
the experimental data for c = 0.1 ppm, Tk = 40 mK, 
s = 1/2 and A = 1.95 ns K 2/3 {A theory = 2.6 ns K 2/3 ). 
This value of c, compatible with the nominal source pu- 
rity, is smaller by two orders of magnitude than the value 
obtained from the fits of energy exchange data on sam- 
ple #1. A similar set of results was also obtained with 
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Cu samples, a material in which the oxyde at the sur- 
face of the films was found to cause dephasing at low 
temperature . Data on energy exchange Eii could 
also be fit with the theory of Goppert et al using 
T K = 300 mK, c = 4.8 ppm, g = 2.3, vJ = 0.1, on top of 
a Coulombic term with intensity K = 0.5 ns _1 meV -1 / 2 
[j^ . This result gives evidence that the anomalous inter- 
actions observed in many Cu wires at B = 0,0] are also 
due to magnetic impurities. Here also, measurements 
of the phase coherence time Q are explained by signifi- 
cantly smaller impurity concentrations (^0.3 ppm). This 
repeated discrepancy on the concentrations deduced from 
the two types of measurements remains an open problem. 
From an experimental point of view, a more quantitative 
test of theory could be obtained in samples with added, 
identified magnetic impurities at a known concentration 

E3- 

To conclude, we have found that anomalous energy 
exchange rates between quasiparticles were strongly re- 
duced by the application of a magnetic field. Moreover, 
the energy and magnetic field dependence of the exchange 
rate can be accurately accounted for by the presence of 
a small concentration of Kondo magnetic impurities |l2| . 
It is worthwhile to compare this result with recent mea- 
surements on Aharonov-Bohm rings, which show that the 
small phase-coherence times found at B = were in- 
creased in a finite magnetic field [28| . All these measure- 
ments indicate that the presence of very dilute magnetic 
impurities is a very plausible candidate to explain both 
extra dephasing and extra energy exchange observed in 
many mesoscopic samples. 
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voret, P. Joyez, C. Urbina and D. Vion. 
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